One-Flavour Hybrid Monte Carlo with Wilson Fermions
نویسنده
چکیده
The Wilson fermion determinant can be written as product of the determinants of two hermitian positive definite matrices. This formulation allows to simulate non-degenerate quark flavors by means of the hybrid Monte Carlo algorithm. A major numerical difficulty is the occurrence of nested inversions. We construct a Uzawa iteration scheme which treats the nested system within one iterative process.
منابع مشابه
Speeding up Lattice QCD simulations with clover-improved Wilson Fermions
We apply a recent proposal to speed up the Hybrid-Monte-Carlo simulation of systems with dynamical fermions to two flavour QCD with cloverimprovement. The basic idea of our proposal is to split the fermion matrix into two factors with a reduced condition number each. In the effective action, for both factors a pseudo-fermion field is introduced. For our smallest quark masses we see a speed-up o...
متن کاملThe Rational Hybrid Monte Carlo Algorithm
The past few years have seen considerable progress in algorithmic development for the generation of gauge fields including the effects of dynamical fermions. The Rational Hybrid Monte Carlo (RHMC) algorithm, where Hybrid Monte Carlo is performed using a rational approximation in place the usual inverse quark matrix kernel is one of these developments. This algorithm has been found to be extreme...
متن کاملAlgorithms for Dynamical Fermions
This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop ‘Perspectives in Lattice QCD’ in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and the...
متن کاملExact 2+1 Flavour RHMC Simulations
Traditionally, “2+1” flavour simulations have been performed using the R algorithm [1], but an exact algorithm is clearly desirable. There also exist exact Polynomial Hybrid Monte Carlo algorithms [2,3], which can be used for 2+1 simulations: however, such algorithms are expensive with regard to memory consumption and for the case of asqtad fermions such an algorithm would be impractical due to...
متن کاملSchwarz-preconditioned HMC algorithm for two-flavour lattice QCD
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32×24, with lattice spacings a ≃ 0.08 fm and at bare current-quark masses as low as 21MeV.
متن کامل